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We present an extension of Adiabatic Quantum Computing (AQC) algorithm for the unstructured 
search to the case when the number of marked items is unknown. The algorithm maintains the 
optimal Grover speedup and includes a small counting subroutine. 

Our other results include a lower bound on the amount of time needed to perform a general 
Hamiltonian-based quantum search, a lower bound on the evolution time needed to perform a 
search that is valid in the presence of control error and a generic upper bound on the minimum 
eigenvalue gap for evolutions. 

■^Z^' . In particular, we demonstrate that quantum speedup for the unstructured search using AQC type 

I ■ algorithms may only be achieved under very rigid control precision requirements. 

o' 

O ■ I. INTRODUCTION AND MAIN RESULTS 

0^ I A. Introduction 

, Quantum computing is believed to possess more computational power than classical computing on certain compu- 
04 tational tasks. For example, factorization of large numbers can be feasible once a quantum computer is built [s^ . The 
basic paradigm which is usually used in the theoretical works on quantum computing is the so called quantum circuit 
model (QCM), see e.g. [13], although the practical realization of it is yet to be found. Farhi and his collaborators 
had proposed the adiabatic quantum computing (AQC) as an alternative, constructive model for implementing 
a quantum computer. It was later realized that from a computational complexity point of view AQC is equivalent to 
all other models for universal quantum computation [1] . 
. Grover's algorithm [23|, originally derived in the framework of QCM, is one of the milestones in quantum computing. 
^ ' The problem it solves can be formulated as following: Given F : {0, 1}" — )■ {0, 1} and knowledge that there exists 
a unique element x such that F{x) = 1, find x. It is clear that classically, one needs to check F{x) for all A^ = 2" 
values of x to find the solution, so the time complexity of doing so is 0{N). Grover's algorithm uses only 0{y/N) 
. , steps to achieve the result, in the framework of QCM. This bound is indeed proven to be optimal for query type of 
. ■ quantum algorithms, see (7[. In the case where there are m (not necessarily 1) values of x for which F{x) = 1, a 
\l [ modified version of Grover's algorithm works in QCM if the number m is known The technique for finding m is 
■ called quantum counting and was developed for QCM in [loj . 

The original motivation of introducing AQC was to derive the physically attainable algorithm for solving optimiza- 
tion problems such as as satisfiability of Boolean formulas by encoding a cost function into the Hamiltonian. One of 
the (very few) models for which AQC had been shown to produce a speedup is Grover's search problem, addressed 
first in [3]. The further works in this direction (e.g. [H, |3l|) considered the original problem of Grover (that is the 
jj] ■ case m = 1). The natural question that arises in this context is whether the Grover type running time, e.g. 0{\/N), 
, is still optimal for a more general class of problem Hamiltonian, characterized by A^ ^ m > 1. That means that 
one wants to derive the analytic lower bound on the runtime of the algorithm, as well as to construct the suitable 
realization of the algorithm for which the running time (being the upper bound on the optimal time) is comparable 
with this lower bound. The partial result in this direction, namely the lower bound for some class of such models, 
was established in (we postpone the more detailed discussion till the next section). 

In this paper, we derive the analytic upper and lower bounds on the amount of time needed to perform a (Hamilto- 
nian based but not necessary adiabatic) unstructured search for the case N ^ m > 1. We also present a lower bound 
on the evolution time needed to perform a search that is valid in the presence of control error and a generic upper 
bound on the minimum eigenvalue gap for the family of the interpolating Hamiltonians used in AQC. In particular, 
we show that in general the 0{\nN/y/N) control precision is necessary in order to achieve quantum speedup over 
classical computation for the small values of m. 
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B. Bounds on the running time in Hamiltonian— based quantum computation 



In the abstract setting of AQC, we are interested in finding the ground state of the given problem Hamiltonian Hp^ 
in the shortest possible time. To this end, we consider a pair of hermitian N x N matrices Hj p, and will assume 
that N ^ 1. Let H{s) be the interpolating Hamiltonian 

H{s):^{l~ f{s))Hi + f{s)HF, (1) 

where / is a monotone function on [0, 1] satisfying /(O) = 0, /(I) = 1. The idea of AQC is to prepare the initial state 
of the system tp{0) in a ground state tpj of the Hamiltonian Hj, and let the system evolve according to the (scaled) 
Schrodinger equation: 

#,(s) = tH{s)Ms), MO) = ^i- (2) 



The adiabatic theorem (AT) of quantum mechanics ensures that under certain conditions (see theorem 11.41 below for 
details) the evolution of the initial state stays close to a ground state of the problem Hamiltonian Hp- For 

AQC to be efficient, the running (i.e. physical) time r in must be much smaller than iV. One then can ask what 
choice of the initial Hamiltonian Hi and the parametrization /(s) minimizes r, and what the optimal value of r is. 

One of the parameters that enters into the upper bound for r in the standard AT is the minimal value g of the 
spectral g ap q (s) between the ground state energy of H[s) and the rest of its spectrum. Consequently, the traditional 
approach [17| to AQC involves the estimation of g. Excluding a very short list of interesting situations for which g 
can be explicitly evaluated (compilation of such examples can be found in (l^'l. it appears to be a hard problem. In 



some instances one can get an idea of what size of g could be by using the first-order perturbation theory [3l[. In 
subsection II CI we present rigorous bounds on the size of the gap for the problem at hand, albeit we don't use them 
explicitly in our study of AQC. 

The main purpose of this work is to obtain the rigorous upper and lower bounds on the optimal running time r for 
a particular class of problem Hamiltonians, satisfying 

Assumption 1. The problem Hamiltonian is of the small rank: 

Rank{Hp) := m < TV . 

This hypothesis is fulfilled in particular for the generalized unstructured search (GUS) problem, see e.g. [lOj. Since 
we are interested in the dynamical evolution of the initial state for which shifting the energy results in the overall 
dynamical phase factor, the above assumption is equivalent to the following condition: Let V denote the largest 
eigenspace of Hp. Then we require that N — dimV <C A^. 

It turns out that for such Hp one can circumvent the standard AT, avoiding the direct estimation of g. We will 
also see that the (nearly) optimal parametrization /(s) is in fact non adiabatic. 

To formulate our results, we need to introduce some notation first: Let ({£^^}^=i) be a set of distinct 

eigenvalues of Hj (respectively Hp), enumerated in the ascending order. It is allowed to the corresponding eigenvalues 
to be degenerate. In what follows, we will denote by Pj (Pp) the eigenprojection of Hj {Hp) onto Ej := El 
{Ep := e(), and by Qj {Qp) the orthogonal projection onto the range of Hj (Hp). To AQC to be meanigful in our 
context we have to impose Ej ^ 0. In the typical setup, Ej — — I. 

Before stating our results, let us note that for AQC to work, it suffices to ensure that '0r(l) has just the non trivial 
overlap with the range of Pp, which we will encode in the requirement ||Pp?/'7-(l)|| > 7 for a "reasonable" 7. Indeed, 
like many quantum algorithms, the AQC algorithm is probabilistic in the sense that it gives the correct answer with 
the probability 7^. The probability of failure can be decreased to the desired value (namely 0{1/N)) by repeating the 
algorithm times. We set 7=1/5 throughout this paper. Another issue that we want to settle is normalization of 
H{s). To that end, we will calibrate Hj p as = \\Hp\\ = 1. Note that without loss of generality we can assume 

that Ep < (since otherwise we can interpolate —Hj and —Hp which only changes the solution tpr of ^ into -ipr). 
We now introduce some parameters in order to formulate our results. Namely, let 5i = WHpijjiW, let 62 — ||PfV'/||; 
and let 6^ = \\Qpipj\\, where Qp is a projection onto Range Hp. Let gp :— E2 — E(. 

Finally we introduce the notion of what we will refer to as a generic Hamiltonian Hj. Given an m-dimensional 
subspace V of C^, the natural question one can ask is what is a distance from the "typical" vector ipj to V. More 
specifically, suppose one has some reasonable probability distribution function for the vectors ^l>I on the unit sphere 
S"^ in (say uniform). Then the expected value of ||</)7|p of the orthogonal projection i^/ of the tpi on V is equal to 
m/N . One can check that the probability of the event {ipi S : |||07|P — m/N\ > am/N} is exponentially small 
in a (see e.g. [isjV Note now that Qpipi is the projection of ipi onto the range oi Qp, which is an m-dimensional 
subspace. We therefore will call Hj generic if its ground state ipi satisfies IIQ.fV'/II = 0{y/m/N). 

Our first assertion is the non-existence result, showing that for any choice of Hj and any function f{s) the running 
time cannot be smaller than r_ defined below. 
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Theorem I.l (The lower bound on the running time). Consider the interpolating family Eq. ([Ij with an arbitrary 
f. Then the running time t_ in Eq. ^ for which \\Pp'ipT-_(l)\\ > 1/5 satisfies 

r > ^^1^, /- ^.<l/5. (3) 

Remark 1. 

1. This result shows that it is impossible to construct the family of the interpolating Hamiltonians H{s) such 
that the evolution of Vl// will have a reasonable overlap with Pp if the running time r is smaller than t_. To 
probe how tight this bound is, one wants to construct a specific family H(s) and the running time t+ for which 
||Pp'0^(l)|| is not small, and make a comparison between t±. We construct such H(s) in the next assertion. 
As we shall see, our bound t_ is not tight {t-/t+ — 1 ^ o{l)), but of the right order of magnitude (meaning 
T-/t^ = 0(1)) in terms of the asymptotic dependence on the small parameter m/N. 



2. For a generic Hj both ^-nd ||PF'0-r|| ^i'g 0{y^ m/N), hence the minimal running time t cannot be smaller 
than 0(v/W 

3. As we will see, the (nearly) optimizing parametrization /(s) is in fact non adiabatic. 



Comparison vi^ith the loannou - Mosca result. 

In [2^ . loannou and Mosca established the lower bound on the running time r for a particular class of problems 
where the initial Hamiltonian Hj is diagonal in the Hadamard basis while the problem Hamiltonian Hp is diagonal 
in the standard basis. Their result is non trivial provided the largest eigenspace of Hj has dimension N — m where 
m <^ VN, and the lower bound they obtained is given by r_ = 0{\fN jra). Since one can always interchange the 
roles of Hi and Hp and shift energy so that the largest eigenspace corresponds to the energy 0, their result can be 
viewed as a slighter weaker version of Theorem II. II for this class of Hamiltonians. 

In the next assertion we construct a specific family i?(s) and determine the runtime t+ for which ||PF'0r(l)|l > 1/5: 

Theorem 1.2 (The upper bound on the running time). There exists an explicit rank one Hi and an explicit function 
f such that ||Pf'0t+(1)|! > 1/5 for 



\EF\\\PpiJi\\' 

and any value C G [1/3,2/3], provided that \\Qp'(l^i\\/ gp = 0(l/lnA^) 
Remark 2. 



1. For a generic choice of Hi one has WQptjjiW = 0{y^m/N), \\PpiJi\\ = 0{^/rrr[/N), where mi = RankPp. It 
implies that t^/t+ = 0(1) for m = 0{1). 

2. This assertion can be viewed as an extension of the result obtained in ^^18] that considered the original Grover's 
search problem in the Hamiltonian-based algorithm. 

3. The interpolating function / in this construction is similar to the one used in Ts'], namely it is a double step 
function. Since / is discontinuous, we prefer to refer to this particular construction as the Hamiltonian-based 
algorithm rather than AQC. 

Note that a-priori the values of Ep and 82 may be unknown. For instance, the value of the overlap 82 has to be 
determined in the GUS problem with the unknown number of marked items. To this end, we prove the following 
auxiliary result: 

Theorem 1.3. Suppose that the value of Ep is knoivn. Then there is a Hamiltonian - based algorithm that determines 
\\Ppipi\\ with l/N"^ accuracy and requires 0((lnA^)^) of the running time. 

Remark 3. 

1. The running time for this sub-algorithm is much shorter than t+, so it does not significantly affect the total 
running time. 

2. A parallel result in the context of the quantum circuit model was establishes earlier in [Toj . 
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C. Gaps in the spectrum of the interpolating Hamiltonian 

Although the size of the gap in the spectrum of H{s) did not play much of the role so far, it is instructive to estimate 
it for the following reason: The size of the gap manifests itself in the adiabatic theorem of quantum mechanics (AT), 
on which AQC is build. The following assertion holds, see e.g. Q: 

Theorem 1.4 (Uniform adiabatic theorem). Suppose that the H{s) is twice differ entiable and hounded family of self 
adjoint operators on the interval [0, 1] that is t -independent, and suppose in addition that 

g := dist(Ai(s),cr(i?(s) \ Ai(s)) > for all s e [0,1] . (5) 

Then the solution ipris) of the IVP (O satisfies 

lim dist {ipr{s), Range Pp) — 0. (6) 

To AQC to be meaningful, one should choose the initial Hamiltonian Hj in such a way that Rank Pj is small. The 
error in the adiabatic evolution (the right hand side of dH])) depends on the size of the gap g, with the rough upper 
bound on the error of the form -S? 1231 . 

Theorem 1.5 (The size of the gap). Let gi := E2 — Ej be a gap between the ground state of the initial Hamiltonian 
Hi and the rest of its spectrum. Let 5^ := \\PiQf\\, where Qp is a projection onto the range of Hp. Then we have 
the following estimate on the size of the gap g in ([5]) ; 

9 < 10<54, (7) 

provided gj > lOS^. 

Remark 4. In fact one can relax the condition gi > 10(^4, but to keep the presentation simple we impose this additional 
constraint. 



D. Robust adiabatic quantum computing 



In this section we propose a necessary technical requirement on the quantum device for AQC Grover's search to 
be successful. The \/]V speedup in AQC algorithm obtained in the tractable problems (c./. [11] for the Grover's 
problem or its rigorous treatment in |23|) relies on a special choice of the parametrization f{s) in ([T]). Namely, it is 
constructed in such a way that f{s) is small at instances {sj} at which the spectral gap g{sj) :— X2{sj) — Ai(sj) of 
H{sj) is small. In the AQC jargon, it is usually referred to as the quantum search by local adiabatic evolution. It is 
interesting to compare this approach with the construction used in Theorem lI.21 where this strategy is pushed to the 
extreme, namely / used there is actually the constant except for the endpoints s = 0, 1 where it jumps. There are 
two practical problems with this approach: 

1. The values {sj} obviously depend on Hp and in particular on Ep (even for the Grover's problem, as the simple 
scaling argument shows). So to choose such an / one has to know the spectral structure of Hp with o{1/Vn) 



precision. This is tacitly assumed in [12|. Note that albeit Theorem II. 2 1 (used in conjunction with Theorem II. 5 1) 



represents an improvement with this regard, it still requires knowledge of Ep. 

2. Even if this technical obstacle can be overcome, the extreme susceptibility of ipri^) to the parametrization / 
poses a radical problem in practical implementation. Indeed, it is presumably extremely difficult to enforce 
/ = for a long stretch of the physical time, as the realistic computing device inevitably fluctuates due to the 
presence of the noise. Some models that try to take into the account the noise were proposed, see e.g. P. [Tl|. 
but to the best of our knowledge all of the existing constructions contain ad hoc parameters and are not derived 
from the first principles. For some interesting rigorous work in this direction that considers de-phasing open 
systems see [6[. 

Another issue [s^ that will motivate our last result below is related to the fact that the adiabatic theorems fall into 
two categories: Those that describe the solutions for all times, including times s G [0, 1], and those that characterize 
the solutions at large times s > 1 where the Hamiltonian is time independent again. Interestingly they give more 
precision for long times. We call the first category, the one that applies to all times, uniform, the second is the long 
time category. 

A representative result from the uniform category is Theorem |L4] above. A characteristic result (see, e.g. (8l.[25l[26j) 
which lies in the long time category is 
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Theorem 1.6 (Long time adiabatic theorem). Suppose that the H{s) is smooth (that is C°° class) and bounded family 
of self adjoint operators with H{s) supported on [0, 1] that is t -independent, and suppose in addition that ([5]) holds 
as well. Then the solution i^ris) of the IVP ([2]) satisfies 

dist (V'r(s), Range Pp) = o(t"") for s>l, (8) 

for any ri G N. 
Remark 5. 

1. In words one can say that starting and finishing the interpolation slowly decreases the error in the adiabatic 
theorem. 

2. There is, in general, no uniformity in n; the term on the right hand side is of order c„t^" where c„ grows rapidly 
with n [c.f. the following discussion). 

3. The distinction between the uniform and the long time AT has an analog in integrals. Suppose that g{s) G 
C°°([0, 1]). Then 



g{t)e''''dt 







o(r-") , if s > 1 ; 
O(r-i) if s e (0,1). 



In the application to AQC it is natural to investigate the dependence of the coefficients c„ in terms of the gap g and 
minimize the run ning time r in such a way that CnT~"' ~ o(l) for some optimally chosen value n. The recent result 
in this direction, [2j, gives r = 0{g~^). For the sketch of the argument that uses (truncated) Nenciu's expansion 
technique [1^1 and leads to the sharper estimate r — 0{g^^\lng\^) see [l3]- One is then tempted to combine the 
starting and finishing slowly strategy with the quantum search by local adiabatic evolution strategy in order to 
minimize the error in the adiabatic theorem. Such analysis was undertaken recently for Grover's search problem in 

mil. 

Inspired by the above discussion, we will assume that in the robust setting for any given moment s inside the 
interval J described below (and which excludes the vicinities of the endpoints s = 0, 1) the value /(s) is greater 
than some small but fixed k > 0. To motivate the definition of J, suppose that the function / lies in the long time 
category, i.e. f is supported in [0, 1] and / is smooth. Let b £ [0, 1] be such that / does not change sign on [6, 1] 
(but it can vanish there). It is not difficult to see that since / is monotone, / has to be concave on [b, 1], hence 
/ < there. Now let us define the interval J for any differentiable function /. Let a — rnin^g[o^i]{/(s) — 1/3}. Let 
b — minsg[o_i]{/(s) is concave on [s, 1]}. We then define the interval J := [a,b] for f € provided a < &, J = if 
a > b, and J := [a, 1] if / ^ C^. To illustrate this notion, consider / S C°°(M) constructed as follows: 



fit) = / gis)ds, gis) = 



, if s ^ [0, 1] ; 
ae^<^ if s e (0, 1) . 

The factor a here is a normalization constant, chosen so that /(I) = 1. We then have 

Rt) = - (1^) ' ^ ^ ' 

so that the only inflection point is t = 1/2. Hence / is convex on [0, 1/2] and is concave on [1/2, 1]. We therefore 
get 6 = 1/2 and J = [/"^ (1/3), 1/2]. The convexity of / on [0, 1/2] implies that f{y) - f{x) > f{x){y - x) for any 
x,y £ [0,1/2]. Choosing x = s, y = 0, we obtain /(s) > f{s)/s for s S (0,1/2]. Since / is monotone, we conclude 
that /(s) > 2/3 on [a, 1/2]. So in this example k = 2/3. 

The utility of the introduction of the interval J is as follows: On the interval [a, 1] the function / is concave, hence 
it satisfies f{y) — f(x) < f{x){y — x) for any x,y £ [6, 1]. In particular, we have 

1 - fit) < f{t){l -t) < fit) for t e [6, 1] , (9) 

the relation we are going to exploit. 

Our last result establishes that in the case of the small rank initial Hamiltonian the robust version of AQC does 
not yield a significant speedup unless n can be made exponentially small: 
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Theorem 1.7 (Robust lower bound on the running time). Suppose that f in Eq. ([T]) is (piecewise) differentiable and 
satisfies f{s)>K>0 for s S J with the interval J defined above. Also, let us assume that Ej = —1. Then, if 
T<Tr^O ( ^^^^i^J , we have 

\{MMm > ^ + 2<5, (10) 

where S — ||Q/Qf||- Hence the running time t for which ||(5fV'i-(1)|| > 1/5 cannot be smaller than Tr- 
Remark 6. 

1. This theorem tehs us that for a generic Hj of the smaU rank the robust running time cannot be smaUer than 
0{kN/ InN). Hence unless the control precision k is on the order of 0{\nN/^/N), AQC is not much better 
than its classical counterpart that solves GUS for r — 0{N). 

2. As we remarked earlier, the requirement Ej — —1 is a very mild one. 



II. PROOF OF THEOREM ITTI 

The proof is based on the following observation: Note that ipi is an approximate eigenvector of H{s) since (-ff (s) — 
(1 — f{s))EiI)ipj — f{s)Hpipj , and the norm of the right hand side is equal to 6i. So by the first order perturbation 
theory, the dynamical evolution of the state ipj given by © will stay close (up to the dynamical phase) to ipi , unless 
the running time t is such that the total variation, given by tSi, is of order 1. The proof below formalizes this 
argument. 

For a solution V'r(s) of ([U, let 

</.,(,,):= e'^WV;,(s), f,{s)=TEi f{r))dr. (11) 

Jo 

Then one can readily check that 4>t{s) satisfies the IVP 

id3r{s) = Tff(s)^,(s), 0,(0) = V/, (12) 

where 

H{s) = {l-f{s)){Hj-EjI) + f{s)HF. 

The factor e^^^'^^'^ is usually referred to as a dynamical phase. 
Let Ur{t, s) be a semigroup generated by H{s), namely 

-idsUr{t,s) = TUr{t,s)H{s); Ur{s,s) = /; t>s. (13) 

Then the solution (pri^) of (IT^ is equal to C/r(l, 0) V"/- On the other hand, 

/-C7,(l,0) = / dsUr{l,s)ds ^ iT [ Ur{t,s)H{s)ds, 
Jo Jo 

hence applying both sides on ipj we obtain 

ipi - (pril) = ir [ Ur{t,s)H{s)iljids . 
Jo 

We infer 

Ui - Mm <r r\\H{s)Mds. 
Jo 

But 

H{s)^i = {{1- f{s)){Hi-EiI) + f{s)HF}i;i = f{s)HFi^i, 



7 



and we get the bound 

Wi^i - 0r(l)|| < t5i f /(s)ds < t5i, 

JO 

where in the last step we used < /(s) < 1. By the triangle inequality, 

- iiPf<^r(i)ii| < iipfV'/ - ppMm < Ui - Mm 

so that 

\\PF4>r{l)\\ < T61+S2. 

On the other hand, by the assumption of the theorem ||Pf-0t-(1)|| > 1/5, hence |lPF0r(l)|l > 1/5. As a result, we can 
bound 

1/5 < tSi + 52, 

and the assertion follows. ■ 

III. PROOF OF THEOREM |L1 

We will choose Hj — —\ipi){ipi\, and a non adiabatic parametrization 

{0 , if s = , 
a^^ifse(0,l), 
1 if s = 1 . 

That means we move extremely quickly (instantly in fact) to the middle of the path, stay there for the time r, and 
then move quickly again to the end of the path. We first observe that regardless of the choice of /(s) in ^ we have 
iPt{s) G Y, where F is a subspace of the Hilbert space, spanned by the vectors in the range of Hp and tpi. Here we 
have used the fact that the range of Hj by the assumption of the theorem coincides with Spanlipi}. Let us choose the 
orthonormal basis {ei}™~![^ for Y as follows: The first m vectors in the basis are the eigenvectors of Hp corresponding 
to {-E/} that differ from zero, and the last vector e„j_|-i is obtained from -0/ using the Gram Schmidt procedure. That 
is, 

Qpipi Qf^^i ^ , ^ 

Gm+l '■= -J— S — JT = , ^„ , = i — Qf ■ 
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We then have 



lem+l-V/ll = ||(5_Fe,„+i - (5_fV/|| + HQFEm+l - Qf->Jji\\ 



2 



= ^3+ (^=f (1-^3) 



< 61+61 (14) 

Our choice of g ensures that 

Here we introduce Pp as the orthogonal projection onto the span of {e^} with e( = e( , and Pm+i the orthogonal 
projection onto Cm+i- Clearly 

p _ QfPiQf 
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We want to compute the matrix elements of the propagator e +Hp) ^^i^ basis {e;}. To this end, we observe 

that in this basis Hp — diag{Ep, . . . , E^, 0), and Hj is a block matrix such that 



EfPi 



S3V* 
S3V Ep 



< iSl\EF\. 



(15) 



where \\V\\ = \Ef\- Indeed, we have 



WQfPiQfW = si, 

1-Pr?j + 1 Pi ^rn+l|| = ^ ^ , 



|P„+i Pi QfW = 63{l-diy^\ 



Together, we obtain that in this basis 

EfPi + Hf - 



D S3V* 
63 V Ef 



< ^SI\Ef\. 



where D = diag{EF, • ■ • , A simple perturbative argument (c/. Duhamel formula (jl7p below) shows that 



with 



K = 



-iarK 



< 3Sl\EF\aT 



(16) 



D S3V* 
S3V Ef 



' D 







■ (53 F*" 







+ 


63V 



To this end, we split K into the diagonal and off diagonal parts: 

K = Ki+ K2 

Let 

then f](s) = —iaTK2{s)D,{s) , and e"*""^^ is given by the following Duhamel formula: 



K2{s) := ^lo^rsK^ J^^^-^^arsK^ 



iyl — iar J K2{s)ds 



Observe now that 



+ i-iary / K2{s)ds / K2{r)n{r)dr 
Jo J{ 



(17) 



(18) 



since e ^2"^^ is diagonal with (1, 1) entry equal to {m + l,m + 1) entry. In fact, [K2{s)]j 



such that Ej — E{. Therefore 



since 



Pf I K2{s)dsP„,+i 
Jo 



PFK2P,n+l = 



\Ep 



= S2\Ef\, 



Pf Pi Pm+i • 



[^2],™+! for all j 



(19) 



To estimate the second term in p^ . we note first that the following bound holds: 
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Lemma III.l. We have 



< 



rgp 9F 



(20) 



This estimate is essentially a content of Lemma 3.3 in [15|. The idea is that, since the spectral supports of Ki 
for Pm+i and QfPf are a distance gp apart, the integral over r has a highly oscillating phase of order rgp- For 
completeness, we prove this lemma below. 

Armed with this estimate and using the fact that i^2(s) is off diagonal, we get 



Pf [ K2{s)ds [ K2{rMr)drP„,+i 
Jo Jo 



= \\PFK2P,n+i\\- P,n+i / ds i^a (r) dr P^+i 

Jo Jo 

<\Ef\52-[ ds f P„^+iK2{r){PF+QFPFMr)P^+idr 
Jo Jo 

< \Ef\52-I^J^ \\Pm+iK2{r)PF\\ dr 

+ max / P„,+iK2{r)QFPF^{r)dr 1, 

where we have used ||f^(r)|| — ||/',„+i|| = 1 and (|T9)) . Applying estimates p9)) and (|20l) . we bound 

Pf I K2{s)ds ( K2{r)n{r)dr P^+i 
Jo Jo 

< \Ef\S2-\\Ef\S2 + ^ + M 
I Tgp gp 



(21) 



Multiplying (fTS)) by Pp from the left and by Pm+i from the right, and using the estimates and (|2T]) . we establish 



FC Pm+l\\ > OLT 



(arf 



Pp [ K2{s)dsP„,+i 
Jo 

Pp f K2{s)ds f K2{rMr)drP,r,+i 
Jo Jo 



aT\EF\62 ■ [1 - ar -i \Ef\S2 
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Tgp gp 

Combining the estimates in (IT4l) . ([16]), and p2|) . the result will follow provided 

8(5^ 261 



(22) 



arlEpl 62 - [l - ar ■ i\EF\52 + — + — 
\ { TgF gp 

> 1/5 + Si + Si + 3Sl\EF\aT. 
Note now that for r = 0(1/62) the above inequality is satisfied for values of r and S3 such that 

aT\EF\S2 ■ (1 - aT\EF\S2) > 2/9, S^/gF = 0(1/ In TV) . 

The result now follows. 
Proof of Lemma I ///._?[ Let 

1 



X := 



2TTi 



Pm+i [Ki - zl)-' K2 (Ki - zl)-' QpPpdz, 



(23) 



where the contour F is a circle {z G C : \z ~ Ef\ = gp /2}. Since 

1 



— ^ (i^i - zl)-'dz ^ Pf + Pm+i , 



10 



one can readily check that 



Hence 



[X,Ki] = Prn+l K2QfPf- 

Pm+iK2ir)QFpFfl{r) dr 



-2i 



Integrating the right hand side by parts, we obtain 



^(^e-^i-^iXe^i'-^i) r!(r) dr. 



P,n+iK2{r) QFpF^ir) dr 



2i 



The first term is bounded in norm by 4-'^, while the second one is bounded by • ||i4r2|l- It follows from 



\X\\ < 



2\\K2\\ 



On the other hand, 11X211 = ^sWVW = \Ef\S3 < S3 by US]), and the resuh follows. 



that 
□ 



IV. PROOF OF THEOREM |Ll 

The algorithm used in the proof is inspired by the mean ergodic theorem and makes use of the fact that the survival 
probability CF{t) = (^A/le**^^ l-i/)/) is directly measurable in AQC framework. We suggest to measure the survival 
probability for a number a times specified below to estimate the overlap 52, and then to count the total running time 
spent on this subroutine. 

Our starting point is a truncated Taylor's expansion for e^: 

fc=0 ^ 



Setting X = pe™, and multiplying both sides by a; ^, we obtain the following relation: 



iwk / r)^\ 

gP(cosu;-l)gjpsinu. ^ g-P ^ ^ Oie~P^\. (24) 



k=0 



If 1 — COS w > g, then the left hand side of Eq. (|24)) is bounded by e~P^ and therefore is smaller than l/N"^, provided 
p = 21nA^/min(l, 5). On the other hand, with such choice of p, the remainder term in Eq. (|24p is bounded by 
0(l/iV^) if L is chosen to be equal to ep. Combining these observations, we get 



f^A^^J i + o(i/iv2), ifc^ = o 

^ kl 1 0{1/N^), ifl-cosa;>(7 ' ^ ' 



k=0 



where p = 2 hiN/ mm{g, 1). 

Now, using the spectral decomposition of Hf, 



N 



where Ei is the «-th distinct eigenvalue of Hf and Pi is the projector onto the spectral subspace associated with Ei. 
Hence 

'^P k 

^""EtT^^^I^'*^''""''"^^^) = {52f + 0{l/N^), (26) 



k\ 



ioT p ~ 2 In N/ min(l, 1 — cos gj?) where we have used Eq. ([25)) . The total running time is J2t=i ^ ~ 0((ln N)'^ ^ 
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V. PROOF OF THEOREM [L5l 

The main tool we are going to use is the so cahed Krein's formula Q for the rank m perturbation of the initial 
Hamiltonian Hj. It gives a characterization of the location of m eigenvalues of the perturbed matrix that differ from 
the spectral values of Hj. Specifically, let A,B be two hermitian matrices, with Rank B = m, and let Q be an 
orthogonal projection onto Range B. Suppose that A is invertible (that is ^ o-(A)). Then Krein's formula tells us 
that 

{A + tB)-^ = (X-i+tBQ)"', (27) 

with 

K := QA-^Q, 

whenever the right hand side of (|27p exists and where K^^ +tBQ is interpreted as acting in the m-dimensional space 
Range B. In other words, is an eigenvalue of A + tB if and only if the m x m matrix K^^ + tBQ contains in its 
spectrum. 

The Krein's formula follows directly from the Schur complement formula, which says that if C is invertible then 



QC-^Q = {QCQ - QCQ (QCQ)-^ QCQ) ' 



where Q := I ~ Q and the inverses on the right hand side are understood as acting on the ranges of Q and Q, 
respectively. 

To apply the Krein's formula in our context, we form a one parameter family 



Ht := Hi + tHp, t= , t e [0, oo) , 

1 — s 



It is then follows that for a fixed t e (0, oo) the eigenvalues of Ht that differ from (t{Hj) are given by the roots of the 
equation 

det {K-\E) + tHpQF) = 0, (28) 

where 

K{E) := Qf {Hi - E)-^ Qp . 

Whenever it is clear from the context that we are working with the operators on Range Qf, we will suppress the Qf 
dependence. 

To analyze ((28)) . we start with the following simple observation: 



Lemma V.l. The matrix K{E) can he decomposed as 

K{E) = k{E) + (29) 

Here the matrix D is positive semi-definite, and is bounded in norm by 1. The matrix K{E) is holomorphic in the 
half plane Re E > Ej — gi/'i, and is positive definite for E G [Ej — gi/2,Ei + <7//2]. Moreover, in this interval we 
have hounds 



^Sl < KiE) < 1; -^-Sl < < 1. (30) 



Proof. We decompose 



K{E) = Qf {Hi - E) ^ Qf 



QfPi {Hi - E) ^Qf + QfPi {Hi - E) ^ Qf . 



The first contribution will correspond to K{E) in ((29| . and the second one to its counterpart in (p9)) . Note now that 
for E e[Ei- gi/2, Ei + gi/2] we have 



3/ 
2 



Pi < Pi{Hi-E) < (2 + |-) Pf , 
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where the the upper bound is a consequence of = 1. Hence we obtain 

^ QfPiQf < QfPi{Hi-E)-^Qf < —QfPiQf- 



4 + 5/ 9i 
Therefore, the first bound in (j30p follow now from 

QfPiQf = Qf - QfPiQf 

and 

< QfPiQf < SIQf- 
To obtain the second bound in ([50)) we note that 

^ {Hj E)-' = {Hi E)-' 

for E ^ a{Hi), and then proceed as above. □ 

In applications to the AQC the parameter 64 is typically extremely small: (5| — 0{l/N). Hence the second 
contribution in (PO)) is small provided \E — Ej\ ^ 64. Therefore for value of E in such intervals, we can first find the 
roots Ei{t) of 



det 



(k'^{E) + tHp) = 0, (31) 



and then estimate \Eiit) — Ei{t)\, where Eiit) are corresponding roots of ([28l) . As we will see, the level crossings or 
the avoided level crossings for Ht occur for values of t such that a pair of eigenvalues Ej.{t), Ei(t) is close to Ej. To 
find these values of t in the first approximation, we fix the value E = Ej in ([5l]l and solve it for t. We have 

Lemma V.2. The equation 

Aet{k-^{Ei) + tHp^ = 0, (32) 

has exactly m+ roots {tj}^^i on (0, 00), where r7i+ is a number of negative eigenvalues of Hp- 
Proof. Let A :— K{Ei), then it follows from previous lemma that < A. Hence 

(A-i + tHp) = (^^-1 ^ A^/^HpA^'^) . 

The right hand side is not invertible for values {tj} such that 

-tj^ e a (a'/^HpA^/^) , 

and the result follows now from Sylvester's law of inertia [2l[. □ 

We are now in position to estimate the size of the gap g from above. Namely, we consider the gaps gj for Ht for 
t = tj. Since t = and H{s) = (1 — s)Ht, we obtain g < y^- < gj. Let 

/ \ 1/2 



4 _ r2 



To get a bound on gj we show that (pS)) has roots in the intervals [Ej — /3,Ei) and {Ej, Ej + f3], a.t t — tj. We then 
infer that gj < 2/3, from which the upper bound in ([7]) follows since (3 < 564. Observe first that by condition of the 
Theorem 01 a-{Hi) Ci [Ei - l3,Ei) = a{Hi) D {Ei,Ei + /3] = (where the latter property follows from the bound 
(3 < 5S4 < gi/2), hence we are in position to use Lemma IV. II We only show that for the first interval, the proof is 
analogous for the second one. 
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To this end, we will denote by sgn(A) the signature of the matrix A. We observe that since in (|29|) is positive 

semidefinite and monotone increasing for the values of E in [Ei — /3,Ei), we have 

sgn{K-\Ei ~Q) + tjHp) < 8gn{k-\Ei) + t^Hp). (33) 

On the other hand, we have 

Ei ^4 



K{Ei) - K{Ei -/?)=/ K'{E) AE > p[ j^^^ - 5l] . (34) 



Ei-I3 



where in the last step we have used (l30l) . Hence 



K{Ei-p) = k{Ei-P) + 



(4 + .g/)2 'V P 



< K{Ej) - 1 ,\,_,, -Sl] I-^] (35) 



< if(i?/) - jl < K{Ei), 

with a choice of (3 as above, and where we have used \\D\\ < 1. We infer 

K-\Ei-f3) + tjHp > K-\Ei) + tjHp, 
and since the matrix K~^(Ej) + tjHp has zero eigenvalue by construction, we obtain 

sgn{K-\Ei - /?) + t.Hp) < sgii{k-\Ei) + t,Hp) . (36) 

Combining and together, we get 

sgn(X-i(£;/ - /?) + iji/F) < sgn{K-\Ei + 0) + iji/F) • (37) 

But the family K^^{E) + tjHp is continuous on [Ej — /3, Ej), hence there should be some value of E in this interval 
for which K^^{E) + tjHp has the eigenvalue 0. ■ 

VI. PROOF OF THEOREM |L7| 

Let us remind the reader that in the context of this assertion Ei ~ ~1. As in the proof of theorem 11.11 

0,(s):= e'/^(^Vr(s), /r(s) = -T I (l-/(r))dr. 

and 

H{s) = (1 - /(s)) (i/, + 1) + f{s)Hp. 
Let us introduce the auxiliary matrix 

B{s) = {f{s)Hp + 1 - f{s) + «:))-i , 

and let 4>{s) — i^i — f{s)HpB(s) tpj, where e is a small parameter to be chosen later. Omitting the s dependence, we 
have 

= -/(I - f)HiHpB - lefHpB . (38) 

That means that away from the m values of s for which B^^{s) has zero eigenvalue, |1^?<?!>|1 is very small, since 
WHiHptpiW < 5"^ . Note now that by fundamental theorem of calculus we have 

(0(1)10.(1)) = imiMO)) + I' -^^{^{s)\Ms))ds, (39) 
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where 0r(s) is defined in But (0(O)|(^r (0)) = 1 and 



I(0(l)l0r(l))| 



(V/|0r(l)> - (H 



-mi)) 



Hp — ei 



< 



Substitution into Eq. ([M)) gives 



Hence Eq. dTO]) will follow if 



1 - liMMm 



< 



n ds 



< 1 



+ 6. 



2Ve 



36. 



To establish the above bound, we note first that 



A. 

ds 



ms)\M''^)) = {Hs)\Ms)) ~ tT{q^{s)\H{s)\Ms)) 



(40) 



(41) 



We bound the first term on the right hand side by ||0|| and the second one by r||i/0||. On the other hand, suppressing 
the s-dependence, we have 

= -fMpB^i - ffB{HF - l)BHFi^i , 



hence 



4 < /I|S|III^F^/Il + 2f\\Bf\\HFibi\ 



where we have used |/| < 1 and / > 0. Let dist(5, z) be an Euclidean distance from the set S to the point z in C, 
and let (j{H) stand for the spectrum of H. Then we can estimate the right hand side further as 



< 



11 + III. 
A. (A,)2 



(42) 



with 



A,(s) distifis)a{HF) , -l + /(s)+ez) 
and where we have used —Qf^ Hf < Qf- Taking the norm from the both sides of we get that 

\\Hcb\\ < WHiQfWW/HfBWWQfM + ^W/HfBWWQfM, 
where we have used HfQf — Hf and ||i//(3F|| < ||Q/Qf||, with the later relation following from 

\\HiQf\\ = WHiQiQfW < IIQ/QfII . 
To estimate \\fHFB\\ we consider three cases: 
1. s G [0, a]: On this interval we can estimate 



(43) 



\\f{s)HFB{s)\\ < max — ^ 

se[0,a] Ae[S) 



< 3. 



(44) 



where we have used a{HF) C [—1, 1] 
2. s e J: In this case, we bound 



(45) 



using theorem's hypothesis. 
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3. s G [b, 1]: Here we estimate 



\\f{s)HFB{s)\\ = \\f{s)HF{f{s)HF + l-f{s) + ei))-^\\ < 1 + 

where in the last step we have used ©. 
Plugging dm) - gSl) into (gSl), we obtain 



< 2 



\\H(t>\\ < (S^+eS) + 
Using and (|T7l) to bound the right hand side of (|1T|) . we get 



In what follows we will use 
Lemma VI. 1. We have bounds 



m 

nAe 



1 r J 

. ^ < -2(TO+l)lne; 
Ac jQ 



^ fds ^ 2(m + l) 



< 3t((5^ +e<5) + 2(m+ 1)5 I -Ine I 1 + — + — ) + - 



(46) 



(47) 



(48) 



(49) 



Integrating both sides of (j4T|) over s and using (j48| and (j49|) . we arrive at 
as 

Hence the required bound in Eq. (|40p follows with the choice e = 10^ (to + 1)6, provided r < —-^f^ where C is some 
generic constant. ■ 



Proof of Lemma [VI. 1[ We derive the first bound, the second bound can be shown analogously. To this end, we observe 
that 



/ 



max 



/ 



< 



E„ea(HF) IfEn + 1 - / - ei| 

/ 



E 



\f{En-l) + l-ei\ 



It follows that 



^ fds fds 

fds 



< (to + 1 ) max / 

'Eel-lA]Jo 



\fE+l-f-ei\ 



where to = Rank Hp. But 



df 



^(f{E-l) + lf + e^ 

^— In [{E - 1)/ + 1 + ^iiE-l)f + ir+e^) 



E - 



1 , / E + y/WT^ 
In 



E-1 \ i + yrr? 
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One can check by taking the derivative that the expression on the right hand side is monotonically decreasing for all 
E € [—1,1]. Since it is also positive and continuous, this term achieves its maximum at E — —1, with the value 




1 



< -2 Ine 



1 + VTT^/ 

for e small enough. Hence the first bound in (|49l) follows. □ 



VII. CONCLUSIONS 



In this work we derived a number of rigorous results concerning Hamiltonian-based quantum search problems that 
satisfy Assumption [TJ Our results include in particular upper and lower bounds on the amount of time needed to 
perform a general Hamiltonian-based quantum search, a lower bound on the evolution time needed to perform a 
search that is valid in the presence of control error and a generic upper bound on the minimum eigenvalue gap for 
evolutions. 

The lower bound on the evolution time is to our knowledge the tightest for AQC type problems. It matches exactly 
the results established in the framework of the quantum circuit algorithm, when applied to Grover's search problem. 

We then construct a specific Hamiltonian-based algorithm with the runtime r = 0{y^N/m) which nearly recovers 
the lower bound, and thus close to be optimal, in the case of Grover's search where the final Hamiltonian is a rank m 
projection. We augment this construction with the Hamiltonian-based quantum counting subroutine, which allows us 
to compute one input parameter 62 to the algorithm. Since the algorithm is robust in a sense that we can allow the 
running time to vary in the large time interval, it is not very sensitive to the ground state energy Ep. As a result, 
just an approximate knowledge of Ep is needed to ensure that the algorithm works. 

While our methods do not hinge on the knowledge of the gap structure of the underlying interpolating Hamiltonian 
H{s), we establish an upper bound on the size of the first spectral gap for this family of matrices. For the final 
Hamiltonian of the Grover's type, i.e. Hp is a rank m projection, the smallest value of the gap cannot exceed 

Lastly, we address the question of the the evolution time for search-Hamiltonians that are also error robust. Namely, 
we obtain the lower bound on the running time when the velocity / is greater than a fixed control parameter n during 
the evolution, excluding the vicinities of the the endpoints s = 0, 1 where it is allowed to be small. We show that the 
necessary control accuracy requirement should be at least OilnN / \/N) for the algorithm to succeed. In particular, 
this result implies that starting slow / finishing slow strategy by itself is not sufficient to make AQC better than the 
classical computer. 

For a general form of Hp for which our Assumption [1] is invalid no general lower bounds on the runtime can be 
obtained, as the examples constructed in fl^ show. It will be very interesting to establish a " typical" lower and upper 
bounds for the random instances of NP-complete problems, discussed in Q. 
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